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ABSTRACT 

The feasibility of a two-way asynchronous (i.e. independently firing) interplanetary laser transponder pair, 
capable of decimeter ranging and subnanosecond time transfer from Earth to a spacecraft anywhere within 
the inner Solar System, is discussed. In the Introduction, we briefly discuss the current state-of-the-art in 
Satellite Laser Ranging (SLR) and Lunar Laser Ranging (LLR), which use single-ended range 
measurements to a passive optical reflector, and the limitations of this approach in ranging beyond the 
Moon to the planets. In Section 2 of this paper, we describe two types of transponders (echo and 
asynchronous), introduce the transponder link equation and the concept of "balanced" transponders, 
describe how range and time can be transferred between terminals, and preview the potential advantages of 
photon counting asynchronous transponders for interplanetary applications. In Section 3, we discuss and 
provide mathematical models for the various sources of noise in an interplanetary transponder link 
including planetary albedo, solar or lunar illumination of the local atmosphere, and laser backscatter off the 
local atmosphere. In Section 4, we introduce the key engineering elements of an interplanetary laser 
transponder and develop an operational scenario for the acquisition and tracking of the opposite terminal. 

In Section 5, we use the theoretical models of the previous sections to perform an Earth-Mars link analysis 
over a full synodic period of 780 days under the simplifying assumption of coaxial, coplanar, circular 
orbits. We demonstrate that, using slightly modified versions of existing space and ground based laser 
systems, an Earth-Mars transponder link is not only feasible but quite robust. We also demonstrate through 
analysis the advantages and feasibility of compact, low output power (<300 mW), photon-counting 
transponders using NASA's developmental SLR2000 satellite laser ranging system as the Earth terminal. 
Section 6 provides a summary of the results and some concluding remarks regarding future applications. 

1.0 INTRODUCTION 

Over 40 international stations, coordinated by the International Laser Ranging Service (ILRS), arc 
dedicated to tracking approximately 25 Earth-orbiting artificial satellites around the clock in support of 
numerous Earth Science applications [Degnan, 1993], These active satellites have ranged in altitude from 
about 300 km to 20,000 km (GPS, GLONASS, Etalon). Apollo 1 1 carried the first of five passive 
retroreflector arrays to the Moon in July, 1969. Because of the great distance between the Earth and Moon 
(-400,000 km) and correspondingly weak signal returns, lunar laser ranging (LLR) has been routinely 
achieved by only three SLR stations over the past 3 1 years. To enhance the weak return signal, all 
successful LLR stations, including the NASA/GSFC funded stations at the McDonald Observatory in 
Texas and on Mt. Haleakala in Hawaii and the French CNES/CERGA Grasse Lunar Observatory in France, 
have been located at astronomical sites with above average atmospheric "seeing" in order to improve 
atmospheric transmission and preserve the narrow laser beam divergence. Furthermore, these stations 
originally employed the largest tracking telescopes and most powerful lasers in the international satellite 
laser ranging (SLR) network. In spite of these advantages, the mean received signal strength per laser fire is 
still well below a single photoclectron because of the familiar dependence of signal strength on range 

for a passive reflector [Degnan, 1993] , In fact, the signal return rates from the Moon are sufficiently low 
that LLR can only be carried out at night and only during certain lunar phases. Near "Full Moon", the solar 
illumination of the lunar surface overwhelms the weak signal return from the target while, at "New Moon", 
no lunar surface features are visible to aid in guiding the extremely narrow laser beam to the target. The 
lunar systems also use special post-detection Poisson filtering techniques which extract the range signal 
from background noise [Abbott ct al, 1973]. 
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In spite of the temporal gaps in lunar range data, the scientific information extracted from centimeter 
accuracy LLR data has been enormous. In addition to precise measurements of the lunar ephemeris and 
librations (deviations from uniform rotation) and of fundamental physical constants, lunar laser ranging 
provides insights into the internal makeup of the Moon and the Earth-Moon interaction and provides 
important tests of General Relativity as well as constraints on its associated metrics [Dickey et al, 1994]. 

The ability to extend these precise range measurements to the inner planets, or to spacecraft in orbit about 
the Sun, would certainly result in a similar array of important scientific results in the fields of planetary 
science, solar physics, fundamental physics, and general relativity. 

Although modem SLR systems are capable of sub-cm ranging accuracies and 50 picosecond time transfer 
between remote clocks over typical near-Earth satellite distances [Degnan, 1993], extending these unique 
capabilities to the inner planets and beyond requires that we adopt a new approach. To overcome the 
prohibitively large R 4 signal loss characteristic of single-ended links to passive reflectors, we must utilize 
two-way laser links, which we will refer to generically as "transponders”. In such systems, which have a 
transmitter and ranging receiver at both ends of the link, the signal strength at each terminal falls off only 
as R 2 . 

The use of microwave or laser transponders had previously been proposed for lunar ranging [Bender et al, 
1990]. Laser transponders on the Moon would make lunar ranging accessible to the smallest of the global 
SLR stations, and the resulting orders-of-magnitudc higher signal strengths would remove the current 
tracking limitations near "New" and "Full" Moon. The NASA Goddard Space Flight Center is presently 
working with colleagues in Japan and Germany to propose a laser transponder on the Selene n Lunar 
Lander Mission [Schreiber et al, 1999], which has been suggested as a possible 2006 follow-on to the 
Japanese Selene lunar lander mission in 2003. 

To date, the most precise ranging to another planet was achieved when the radio telescopes of the 
NASA/JPL Deep Space Network (DSN) were in communication with the Viking Lander on Mars. This 
two-way microwave link typically resulted in several meter precision measurements, with 3 meters 
precision being the best internal consistency ever reported over a short tracking interval [Lemoine, 2000]. 
The absolute accuracy of these interplanetary measurements, however, is limited by the presence of the 
interplanetary solar plasma, which, like the Earth's similarly charged ionosphere, can have a significant, 
and largely unknown, impact on the propagation delay at microwave frequencies. Optical frequencies, on 
the other hand, are much too high to be affected by the charged solar plasma. Furthermore, the shorter 
optical wavelengths characteristic of lasers allow the transmitted energy to be propagated in highly 
collimated (few arcsecond divergence) beams and therefore require much smaller collecting and 
transmitting apertures at both ends of the link Of course, mutual high accuracy pointing between the space 
and ground terminals is required to take advantage of the narrow laser beam divergence. 

The feasibility of interplanetary laser transponders was introduced in the Proceedings of the last two 
International Workshops on Laser Ranging [Degnan, 1996; Degnan et al, 1998]. The present paper seeks to 
review this earlier material, develop some additional models, and provide a detailed analysis of an Earth- 
Mars link In Section 2 of this paper, we will describe two types of transponders (echo and asynchronous), 
introduce the transponder link equation, describe how range and time can be transferred between terminals, 
and discuss the potential advantages of photon counting asynchronous transponders for interplanetary 
applications. In Section 3, we will discuss and provide mathematical models for the various sources of 
noise in an interplanetary transponder link including planetary albedo, solar or lunar illumination of the 
local atmosphere, and laser backscatter off the local atmosphere. In Section 4, we discuss the key elements 
of an interplanetary laser transponder and develop an operational scenario for the acquisition and tracking 
of the opposite terminal. In Section 5, we use the theoretical models of the previous sections to perform two 
Earth-Mars link analyses over a full synodic period of 780 days. We demonstrate that, using slightly 
modified versions of existing space and ground based laser systems, an Earth-Mars transponder link is not 
only feasible but quite robust. We also demonstrate through analysis the advantages and feasibility of 
compact photon-counting transponders using NASA’s developmental SLR2000 satellite laser ranging 
system as the Earth terminal. The latter system uses few mm long, passively Q-switched Nd:YAG 
microchip lasers to generate subnanosecond pulses [Degnan and Zayhowski, 1998]. Section 6 provides a 
summary of the results and some concluding remarks regarding future applications. 
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2.0 BASIC TRANSPONDER THEORY 


2.1 Transponder Types 

Transponders fall into two basic classes, i.e. echo and asynchronous. The type typically used in underwater 
acoustic sounders or microwave aircraft transponders is the echo transponder in which a pulse emitted 
from terminal “A” is detected by terminal ct B” which then generates a response pulse subsequently detected 
by “A”. The delay between the received and transmitted pulse at the secondary terminal, "B", is usually 
either known a priori through careful calibration or controlled via active electronics and is subtracted from 
the observed roundtrip time before computing the target range. Alternatively, the delay can be measured 
locally at terminal "B” and transmitted to the primary terminal , "A", via a communications link. The signal 
return rate at the primary station is then equal to the fire rate of the laser multiplied by the joint probability 
that pulses are detected at both ends of the link This simple echo approach works very well when the 
roundtrip time -of- flight is relatively short and there is a high probability of detection at both ends of the 
link, i.e. when both the uplink and downlink signal is reasonably strong and pointing uncertainties are small 
relative to the transmitter divergence. However, in situations where the light transit time is relatively long 
and the probability of detection is small at one or both ends of the link, it is worthwhile considering an 
alternative approach, i.e., the asynchronous laser transponder . 

In an asynchronous transponder, the two terminals independently fire pulses at each other at a known laser 
fire rate,^ as illustrated by the timing diagram in Figure 1. For an Earth-Mars link, for example, the Earth 
terminal records the times of departure of its own transmitted pulses ( fe; ) as well as the times of arrival of 
pulses from Mars (fo)and vice versa. In a high SNR system, the pulses arrive at rougjily the laser fire rate 
whereas, in low SNR systems, the pulses may arrive intermittently. The departure and arrival times 
measured at each terminal are then communicated to, and properly paired at, an Earth-based processor 
which then calculates a range and clock offset between the two terminals for each set of two way 
measurements occurring within a reasonably short time interval. In a manner analogous to the familiar 
lock4n amplifier , a priori knowledge of the repetition rate at the opposite terminal, as well as the range and 
range-rate between terminals, allows the special correlation range receiver (CRR) to extract the relatively 
scarce, but temporally coherent, signal photons from the more prevalent, but randomly distributed, 
background noise photons [Degnan. 2000a]. 

2.2 TWO-WAY MEASUREMENT OF RANGE AND CLOCK OFFSET 

If we assume that both lasers are firing at a common rate there are four possible measurement outcomes 

within a given "transponder cycle" or period. If Pq A and Pq B are the Poisson probabilities that a pulse is 
independently detected at terminals A and B respectively, each outcome has its own probability of 
occurrence: 

(1) No incoming pulses are detected at either terminal (null event) = (1-P D A )(1-P D B ) 

(2) An incoming pulse is detected only at Terminal A - P d a (1-Pd B ) 

(3) An incoming pulse is detected only at Terminal B = P d b (1-Pd A ) 

(4) Incoming pulses are detected at both ends of the link = P D A Pd 

Outcome (1) is a null event ; (2) and (3) are one-way events originating from opposite terminals within the 
same transponder cycle respectively; and (4) is a two way event where range pulses are observed at both 
ends of the link within the same transponder cycle. Any two one-way events originating from opposite 
terminals, whether they occur in the same transponder cycle or not, allow the post facto calculation of an 
instantaneous spacecraft range and spacecraft clock offset via the equations: 

R = —{ ( BA = “*4l) + (*S2 -O] W 
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where the intervals (t Al - t AX ) and (t B2 - t BX ) are measured by Terminals A and B respectively. In (1) 

and (2), R and t are the instantaneous range and clock offset at the point in time when the “photon world 
lines” marked and Iba in Figure 1 cross each other as can be deduced by simple symmetry arguments. In 

• 

(2), the small correction term, R /c, corresponds to the instantaneous range rate between the Earth station 
and the spacecraft divided by the speed of light and can be estimated from a priori planetary ephemeredes, 
the doppler shift in the microwave communications link, or iteratively solved for from the laser range data. 
Recent analyses for the Earth-Mars [Lemoine, 2000] and Earth- Vesta [Giorgini, 2000] links using 
ephemerides provided by JPL have both shown that, while the error in a priori interplanetary range may be 
as much as several tens of kilometers, the error in the range rate is quite small, typically a few mm/sec. 

As will be discussed in Section 4, frequent one-way events provide a means of acquiring and tracking the 
opposite terminal in both angle and time (range), thereby improving the probability of a near simultaneous 
two-way event. The ability to use any two one-way events in the same or separate transponder cycles to 
compute range and clock offset marks the principal technical advantage of asynchronous over echo 
transponders. In low SNR interplanetary experiments, where the probability of detection can be small on 
both ends of the link, the number of one-way events will far outnumber the two-way events required by an 
echo transponder. This will be illustrated in Section 5 for a simulated Earth-Mars link. 

2.3 ACCURACY AND VARIANCE OF THE RANGE AND CLOCK OFFSET MEASUREMENTS 


Since pulse transit times are on the order of several minutes, spacecraft clock frequency offsets and 
instabilities dominate the systematic and random interplanetary ranging errors respectively whereas, for 
conventional satellite laser ranging, the light roundtrip transit time is typically measured in milliseconds, 
and uncertainties in the atmospheric propagation delay dominate the systematic error at the sub-cm level. 
From Equations (1) and (2), the systematic errors in the measured range and clock offsets due to a constant 
frequency offset error in the Earth and/or spacecraft clocks from some nominal frequency f c are given by 
the equations 
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Wc see from the latter equations that the fractional error in range is equal to the average of the fractional 
frequency offsets in the two clocks whereas the fractional error in the clock offset determination is 
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proportional to the difference between the fractional frequency offsets in the two clocks. The variance in 
the range and offset measurements is in turn given by 
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for two independent free-running clocks. 

Thus, the precision and accuracy of the single shot ranging and time transfer will be limited by the less 
accurate of the ground and spacebome clocks. If the spacebome clock had the stability of a good rubidium 
standard (5 x 10' 13 ), sub-decimeter level single shot range uncertainties would result for distances on the 
order of 1 AU (-150 Million Km). With respect to single shot time transfer, the latter clock would 
introduce a timing jitter on the order of 250 picoseconds (< 4cm). However, if both clocks were of maser 
q uality (~1 x 10 13 over time intervals of several minutes), clock errors would introduce submillimeter 
errors over distances on the order of 1 AU and range accuracy would then be limited, as in conventional 
SLR to artifical satellites, at the sub-cm level by uncertainties in the atmospheric propagation paths. 

One potential enhancement that has yet to be explored thoroughly is the degree to which an inferior 
spacecraft clock can be "disciplined" by an Earth-based maser during transponder operations via the laser 
time transfer process. This would be somewhat analogous to the present-day "disciplining" of quartz crystal 
and rubidium frequency references via frequent comparisons with the 1 pps output of GPS timing receivers. 
The latter receive timing signals from a maser-controlled satellite constellation and results in a long term 
frequency stabilization and performance of the slaved clock on a par with the master clock. 

2.4 TRANSPONDER LINK EQUATION 

The mean number of signal photoelectrons recorded by the terminal B receiver on a single laser fire is 
given by the transponder link equation 


b C ab E a A b (7) 

R 2 

where Ea is the pulse energy transmitted from terminal “A” , As is the area of the receiving telescope at the 
receiving terminal “B”, R is the distance between the two terminals, Cab is a transponder constant given 
by 
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rj q B is the detector B quantum efficiency, h v is the laser photon energy, T A and T B are the one-way zenith 
atmospheric transmissions at Terminals “A” and “B” (appropriate for two terminals on planetary surfaces), 
9 a and 9b are the local zenith angles for the opposite terminal, Q? is the transmitter solid angle for laser A, 
and 7 ] r B is the optical throughput efficiency of the receiver B optics respectively. The order of the subscripts 
in the transponder constant determines the direction of the link, i.e. Cab is the link from terminal A to B. 
The mean signal at the opposite terminal is obtained by simply interchanging A and B in the above 
equations. 

2.5 "balanced" transponder links 

We can choose to write (7) in terms of the average laser power, i.e. 


n . 
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where f qs is the Q-switching frequency (laser repetition rate) in Hz and we refer to PaA b as the Mixed 
Power-Aperture Product to reflect the fact that it contains properties of two different terminals. Although 
not essential, the use of a common laser fire rate in the double-ended transponder link conceptually 
simplifies the instrument. From (9), the ratio of the mean signal strengths at the two terminals is given by 
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where, from (8), the approximation holds if the laser beam divergence and detector and optical throughput 
efficiencies are roughly equal at both ends of the link ( C AB = C ^ ) . Furthermore, the above ratio is 
approximately unity if the Power-to-Aperture Ratio is the same at both ends of the transponder link; i.e. 
P a /A a = Pb/Ab. We will refer to this as a balanced system. Note that the signal detection rate is 
proportional to the Mixed Power- Aperture Product, PaAb, and that, for a balanced system, the frequency of 

measurement at both terminals is approximately equal, i.e. . Thus, if it becomes necessary to 

conserve limited spacecraft resources while still maintaining a balanced link with a fixed sampling rate, one 
can increase the power and telescope aperture at the Earth station proportionally to maintain the same 
Power-to- Aperture ratio while simultaneously reducing the transponder laser power and receive aperture at 
the remote terminal by the same factor. 

2.6 MAXIMIZING THE SAMPLING RATE THROUGH SINGLE PHOTON DETECTION 

The mean number of range returns per second recorded at either terminal is given by the product of the 
laser fire rate and the probability of detection for that terminal, i.e. 
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where we have solved for in (9). The probability of detection in the low signal limit is given by Poisson 
statistics and depends on both the mean signal strength, n,, and the detection threshold, n t . Thus, for 
terminal “B”, the maximum range sampling rate is 
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In Figure 2, we plot the normalized detection rate, fjfmax and demonstrate that the parameter/*** 
corresponds to the maximum detection rate in the limit of small mean signal strength and single 
photoelectron detection. We draw the following conclusions from the figure: 

For a given Mixed Power-Aperture Product, the maximum range sample rate is obtained by using 
single photon detection and a high repetition rate laser with low energy per pulse such that the mean 
signal count received per laser fire is ii* < 0. 1 . 

For higher detection thresholds, the frequency of signal returns peaks at a value lower than / m4X for 

mean signal counts near the threshold value (e g. at rt s = 1.8 pe for n t = 2 pe and at n* = 3.3 pe for n t = 

3 pe, etc ). To achieve comparable return rates at thresholds higher than one photoelectron, the 
instrument power-aperture product would have to be increased significantly (e.g. by factors of 3.4 and 
5. 1 for thresholds of 2 and 3 pe respectively). 

In addition to the sampling rate advantage, there is much to be gained from a size, weight, and power 
efficiency standpoint if one operates in single photon detection mode at both ends of the link. This is 
especially relevant when one considers interplanetary distances of several hundred million km and the R 
dependence on signal strength. Severe size, weight, and prime power limitations on interplanetary 
spacecraft typically preclude either the use of large and powerful laser transmitters to enhance the signal 
strength at the Earth station or large telescopes to increase the signal strength at the spacecraft. Fortunately, 
as mentioned previously, one can place as much of the transponder link burden on the Earth station as 
necessary (within reasonable engineering or economic limits) in order to achieve a "balanced” and robust 
link While Figure 2 provides us with some interesting food for thought, we must have a viable approach to 
extracting single photon signal returns from the relevant noise background. 

3.0 NOISE SOURCES AND MODELS 

Noise count rates can vary widely depending on whether the transponder is operating (1) in cruise phase 
between planets, (2) in orbit about the target planet, or (3) from the surface of a planet containing a 
scattering atmosphere. In this and subsequent sections, we will use the term "planet” in its broadest sense 
to include not only the major planets but also the minor planets (asteroids), moons of major planets, etc. 

The following three noise sources are present in all of the aforementioned operational scenarios: detector 
dark counts, solar radiation scattered from the surface and atmosphere of the planet being viewed by the 
receiver (planetary albedo), and the residual stellar background in the receiver field of view (FOV). The 
number of counts from the star background within a narrow receiver FOV on the order of 100 prad (20 
arcseconds) are expected to be relatively small compared to those induced by planetaiy reflections or 
emissions. Mars under maximum solar illumination, for example, appears as a magnitude -2.5 star when 
viewed from Earth [Zissis, 1993]. For a planetary lander operating in local daylight, we must add the noise 
background caused by solar scattering in the local atmosphere. During local night operations, there may be 
lesser background rates due to the irradiance of the local atmosphere by nearby planetary moons. Back- 
scattered laser radiation produced by the outgoing pulse is an additional time-dependent source of noise and 
is applicable to both day and night operations from the planetary surface. 

At the asteroid belt on the outer limits of the inner Solar System, the Sun is never more than about 25 
degrees away from the transponder line-of-sight to Earth. Thus, direct solar illumination of the transponder 
optics and the resulting scatter within the instrument is another potential source of background noise. 
System baffling and stray light rejection are therefore important considerations for any practical instrument, 
but the noise background count rate can be reduced to acceptable levels through the use of solar shields, 
careful spectral and spatial filtering, and aggressive stray light control. 

Compared to other sources of noise to be discussed in subsequent subsections, dark count rates in the 
visible detectors typically used in laser ranging tend to be relatively low (10" to \0 A counts/sec). 
MicroChannel Plate PhotoMultiplier Tubes (MCP/PMT’s) typically have far lower dark count rates but 
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their quantum efficiencies are often significantly reduced relative to Silicon-based Avalanche Photo-Diodes 
(APD’s). 

3.1 SOLAR AND LUNAR SCATTERING IN THE LOCAL ATMOSPHERE 

Since the effective thickness of the atmosphere, h B (-20 Km for Earth) is small relative to the volumetric 
mean radius of the planet, R B (6378 Km for Earth), we can accurately represent the atmosphere (except 
very near the horizon) by a stratified rectangular slab with the station centered on the bottom face. In our 
model, the atmospheric extinction coefficient, (j. (z), will be assumed to be a function of the vertical height 
only, i.e. no horizontal gradients are assumed. The atmospheric transmission at zenith is then given by 


Tg = Tg{h B s ,<x>) = exp 



(13) 


where h s B is the height of the B terminal station above some reference height (eg. sea level on Earth). 

The solar irradiance can be viewed as a plane wave incident on the top plane of the atmosphere at a solar 
zenith angle, ft, as viewed from the terminal on Planet B. If we ignore the slight bending of the solar rays 
as they propagate through the vertically stratified atmosphere, the transmitted irradiance to a differential 
volume at an arbitrary height above sea level, z, is given by 


Affexp -sec 0 s jdz' ju B (z') = N B [T B (z,cc)Y" d: 


(14) 


where N* is the exoatmospheric solar irradiance at Planet B. At the 532 nm wavelength, the solar 

irradiance at Earth has a value 0.2 watts/m 2 -A° and falls off as MRb where R B is the distance of Planet B 
from the sun expressed in Astronomical Units (1 AU — 150 million km). 


The receiver has a solid angle field of view (FOV) given by 


Of = nOlg 


(15) 


where 0 r e is the half angle of the far field cone viewed by the receiver. When the opposite terminal is at a 
zenith angle 0a the scattered solar power collected locally by a receiver with collecting aperture Ab and 
spectral bandwidth AA, B is given by the following integral 
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where 

s = (z - h B s ) sec 0 B (16b) 

is the range to the differential volume, is the volumetric scattering coefficient at a height z above sea 
level, and we have assumed isotropic scattering. For wavelengths far from absorption lines, the atmospheric 
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extinction is due entirely to scattering so we can set j3 B (z) = n B (z). Using (13) and (16), the rate at which 
solar background light is converted into photoelectrons by receiver B is 
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where q q B is the detector quantum efficiency, q r is the receiver optical throughput efficiency, and hv is the 
photon energy at the laser wavelength. In order to evaluate the final integral in (17), we define the new 
variable 
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Substituting (18) into (17) and evaluating the trivial integral yields 
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(19) 


It should be pointed out that the atmospheric noise model makes no assumptions regarding the distribution 
of scatterers with altitude. The model does assume no horizontal gradients, however, and therefore depends 
only on the zenith transmission, T B , between the station and the "top" of the planetary atmosphere and the 
local zenith angles of the Sun, ft, and the opposite terminal, 6b. Note also that as 7g — >1 (no scattering), the 
background count rate due to solar scatter off the atmosphere correctly goes to zero for all values of Os and 

ft. 

When the solar and spacecraft zenith angles are equal (ft = ft), the bracketed [] term in (19) becomes 
indeterminate, but we can use L’ Hospital’s rule to obtain 
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for the special case when co = sec0 5 - secOs — 0. This approximation also holds in the limit of a highly 
transmissive atmosphere for arbitrary Os and Ob yielding an approximate expression 


”ls = 


vX N’(M B )A B Cl‘ 


hv 


An 


{t™°° In t ?°- } 


( 21 ) 


9 



which is independent of the solar zenith angle. In the case of a planetary lander communicating with an 
Earth station, (19) can clearly be used to assess noise at the opposite terminal by interchanging the 
superscripts and subscripts, A and B. 

Plots of (19) and (21) show the following general results: 

• For all values of the spacecraft zenith angle, Ob, or zenith transmission, T B , the noise rate due to local 
atmospheric scattering of solar radiation is always greatest when the Sun is directly overhead {0 S — 0) 
and decreases monotonically to zero as the Sun drops below the horizon. 

• For higher values of the zenith transmission and all solar zenith angles, the noise count rate grows from 
a minimum at zenith and reaches a maximum at very high spacecraft zenith angles {Ob >80°). This is 
illustrated in Figure 3 where we have chosen T B — 0.7. For very high atmospheric transmissions 
approaching unity, the noise count rate at most spacecraft zenith angles approaches zero because of the 
low density of scatterers contained within the viewed atmospheric volume except very near the horizon 
where the long atmospheric path length and correspondingly larger viewed volume compensates for 
the low density of scatterers. 

• For low values of the zenith transmission (not shown in Figure 3) and all solar zenith angles, the noise 
count rate drops off monotonically from a peak at zenith. This is due to the poor penetration of solar 
radiation at high zenith angles. 

• The approximate expression in (2 1 ) works best at spacecraft zenith angles below 60 degrees. 

At full Moon, the lunar spectral irradiance outside the atmosphere has the value 4.8 x 10' 7 W/m 2 -A° 
at a wavelength of 532 nm [Zissis, 1993], or over six orders of magnitude less than the solar value, which 
can be substituted into the preceding equations to determine the worst case reflected lunar light scattering 
from the atmosphere. 

3.2 PLANETARY ALBEDO 

We can use the results and techniques of the previous subsection to develop a model for the photoelectron 
count rate resulting from solar scattering from the viewed planet (planetary albedo). In Section 4, we 
propose to use these photoelectrons in the initial acquisition of the opposite terminal by imaging the host 
planet onto a pixellated imaging detector such as a CCD array. In the ranging portion of the receiver, these 
count rates also represent an inescapable source of background noise, which the laser signal from the 
opposite terminal must overcome. 

The background rate from planetary albedo can consist of two components - surface and atmospheric (if 
any) scattering. The worst case background occurs when the two communicating planets and the Sun fall 
on or near a linear configuration, which also coincides with the points of minimum and maximum 
interplanetary range (see Section 5). Near the point or closest approach, the terminal on the innermost 
planet is in darkness and ranging to the terminal on the outer planet which is bathed in daylight and cannot 
see the inner planet against the blinding background of the Sun. Near maximum range, the two terminals 
are both in daylight and the transponder link is blocked by the Sun itself. Thus, transponder operations can 
only occur when the aspect angles, a A and a B (see Figure 4a), are sufficiently far from the values 0° or 180° 
to avoid direct observation of the Sun by one or both terminals. The fractional disk area illuminated by the 
Sun drops monotonically from 1 to 0 as the aspect angle varies between 0 and 180 degrees as shown in 
Figure 4b. The fractional disk area viewed by transponders on Earth and Mars are plotted as a function of 
synodic phase in Figure 4c (see Section 5). 

In modelling planetary albedo contributions to transponder noise, we assume a spherical planet surrounded 
by a relatively thin atmospheric shell. The noise count rate in receiver B due to solar radiation scattered 
from the surface of Planet A is given by the equation 
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where we have assumed Lambertian scattering from the surface and included the effects of the atmospheric 
transmission at both planets but ignored the slight geometric bending of the solar rays in the atmosphere. 

In (22a), p A (0, <p) is the reflectance of Planet A’s surface at the operating wavelength, r A is the volumetric 
mean radius of Planet A, 0 B is the local zenith angle of Terminal A as viewed from Terminal B, 


sec y - [sin 9 sin a B cos <p + cos 9 cos a B ] ( 22b ) 

and ob is the aspect angle , defined as the angle subtended by two vectors, both originating from Planet B, 
directed toward the Sun and Planet A respectively. The two-dimensional integral in (22a) is over the solar- 
illuminated surface area within the receiver FOV. If Planet A is sufficiently far away so that it is contained 
totally within the Terminal B FOV, the integral is over the area where the hemispheres in the direction of 
the two aforementioned vectors overlap. If the opposite planet is contained entirely within the receiver 
FOV and the atmosphere is either highly transmissive or non-existent (T A = 1), (22a) can be approximated 
by 
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has a value between 0 (a A = n) and 1 (a A = 0), and, since the planets spin on their axis, we have substituted a 
globally-averaged value for the surface reflectance, p A . 

Following an analysis similar to that in Section 3. 1, the noise count rate in receiver B due to solar power 
scattered from the atmosphere of Planet A can be shown to be 
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where h A is the height of the atmosphere above Planet As surface. In (24), we have used the spherical 
coordinate versions of (13) and (14) and again assumed that there is little or no atmospheric absorption at 
the laser wavelength and that the volumetric scattering coefficient, /?, can be replaced by the volumetric 
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extinction coefficient, /u. In our thin atmosphere approximation, we have h A « r A and can therefore 
approximate fby r A in the last integral in (24). Using (18), we can then evaluate the radial integral to 
obtain 
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which, in the absence of a scattering atmosphere on Planet A (T A =1), correctly gives zero for the 
atmospheric contribution to the background For a highly transmissive atmosphere, we can approximate the 
exponential in (25) by the first two terms in a Taylor expansion, which yields the approximate expression 
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where the geometric factor for atmospheric scattering, 
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again takes on values between 0(04 = 71 ) and 1 ( a A -~ 0 ). 

Summing (23a) and (26a) gives an approximate expression for the noise count rate for the total planetary 
albedo, i.e. 
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where the subscript sa refers to surface and atmosphere. 


3.3 LASER BACKSCATTER FROM THE LOCAL ATMOSPHERE 

Unlike the sources of noise discussed previously, which are all constant in time, the contribution of laser 
backscattering to the overall noise count falls rapidly with the time after pulse emission, x. At large 
spacecraft zenith angles (0 B ~ 70°), the effects of laser backscatter at the Earth station can be significant for 
time scales up to 100 psec, which corresponds roughly to the time it takes laser photons to reach the 
uppermost layer of the atmosphere and return at low elevation angles. Thus, photons arriving from the 
opposite terminal during this early segment of the local laser fire interval will be competing with an 
additional time-dependent noise source. The effects of laser backscattering can often be suppressed in 
ranging systems by employing a bistatic design , i.e. one having truly separate transmit and receive optics 
with parallel but not coaxial paths. Where common telescope optics are used for the transmit and receive 
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beams, an alternative approach is to spatially isolate the transmit and receive beams through a combination 
of “aperture sharing” and aggressive stray light suppression. Bistatic receivers have the additional virtue 
that "transmitter point-ahead" can be more easily accommodated (see Section 4.2). 

For the sake of completeness, the backscattered laser power collected by a coaxial, monostatic receiver B is 
given by the lidar equation 
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where s is again the range to the scattering volume, c is the speed of light, and E B is the transmitted laser 
energy. For the stratified and non-absorbing atmosphere considered previously, the rate of back-scattered 
photoelectron generation at the detector is equal to 
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Experimentally, one finds the following height dependence for the extinction coefficient [RCA, 1968] 


Mb 09 = Ml exp (-Z /hf c ) (30) 

where fi B ° is the extinction coefficient at the nominal lowest-lying surface (sea level on Earth) and hj* is 
the atmospheric scale height for Planet B. Substituting the latter expression into (13) yields the relation 
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Using (30) and (3 1) in (29), we can write 
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We can now choose to express (32) as a function of the time from laser fire, i.e. 


T - ■ 


2s 2 sec 0 \ 


c c 

to obtain our final result 
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which, due to the logarithmic term, also goes to zero in the absence of a scattering atmosphere. 
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4. ACQUISITION OF THE EARTH TERMINAL 


A block diagram of a dual mode Microlaser Altimeter and Transponder (MAT) is provided in Figure 5. 
When in orbit about the planet, the instrument ranges to the surface and provides topographic maps 
[Degnan, 2000b], The same instrument can function as a transponder at any point within the mission, i.e. 
during interplanetary cruise phase, in orbit or on the surface of the planet. In transponder mode, acquisition 
of the Earth terminal requires an initial search within a three-dimensional volume bounded by the initial 
an gular pointing uncertainty and the uncertainty in the time of arrival of the pulse from the opposite 
terminal. As mentioned previously, the latter determines the width of the range gate and has two 
components - the uncertainty in the a priori knowledge of range between the two terminals as derived from 
ephemerides and the uncertainty in the laser fire time at the opposite terminal. As will be described in the 
following subsections, acquisition of the opposite terminal is most easily accomplished in two steps - first 
in 2-D angular space and then in 1-D range space. 


4.1 INITIAL ACQUISITION OF THE EARTH STATION IN ANGULAR SPACE 

The angular search for the Earth terminal is aided by a sensitive CCD array capable of imaging the Earth, 
Moon, and nearby stars within a nominal 1° x 1° degree field of angular uncertainty as in Figure 6a. For a 
transponder mounted to the spacecraft body, this level of angular uncertainty is representative of the error 
associated with the pointing of a spacecraft from orbit or during interplanetary cruise phase. It is also 
inclusive of the expected angular error in the case where the transponder is assumed mounted to a meter- 
class K-band microwave communications dish communicating with Earth. Using pointing corrections to an 
independent two-axis transponder gimbal mount of limited angular range (<2°), or alternatively a pair of 
Image Motion Compensators (IMC's) [McElroy et al, 1977], to implement the fine pointing of the receiver, 
the system computer can center and hold the Earth image in the CCD array. Space-qualified, high 
sensitivity CCD cameras with up to 2048 x 2048 pixel resolution are readily available and yield a 8.8 prad 
single pixel resolution for a nominal l°xl 0 array FOV. Since the full Earth disk subtends an angular width 
between 34 and 163 prad (i.e. 4 to 19 pixels across) from Mars at its farthest and closest points from Earth 
respectively, the center of the Earth image can be well resolved at the sub-arcsecond level. In planetary 
orbit or during cruise phase, the transponder CCD sees a rather bright sunlit Earth against a dark 
background except when ~ 0° or 180°. However, even when the Earth’s “nightside is largely directed 
toward the transponder, there can be sufficient forward scattering of solar light by the Earth’s atmospheric 
rim for detection [Mallama, 1998], 

In order to further quantify our ability to detect the Earth under different operational scenarios, such as a 
lander mission, let us assume that both terminals contain a square CCD array containing pixels 
Assuming that the Earth image does not move substantially during the CCD integration time, the number of 
pixels illuminated by the partially sunlit Earth image (Planet A) is approximately 
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where Qccd is the solid angle viewed by the Terminal B CCD array, /a represents the fraction of the Earth 
disk area illuminated by the Sun where 
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and Oa is once again the aspect angle introduced in Section 3.2. The number of illuminated pixels as a 
function of the synodic phase is illustrated in Figure 6c for terminals operating at Earth and Mars. The 
Moon has the same aspecrt angle as Earth when viewed from Mars, but the image occupies fewer pixels 
due to its smaller size. 

The worst case noise scenario is for a transponder on the planet's surface operating in daylight in the 
presence of a scattering atmosphere. The resulting background count rate is distributed over all of the pixels 
whereas the solar radiation scattered from the opposite planet is distributed only over N A pixels. We can 
now use (19) and (27) to obtain an expression for the CCD contrast of the planetary image at terminal B 
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where the approximation holds for either a highly transmissive atmosphere or when the spacecraft and 
zenith angles are comparable. Figure 6d plots the contrast of the Mars, Earth, and Moon images against the 
scattering of solar radiation by the local atmosphere during daylight operations as a function of the synodic 
phase. 

Once the Earth is centered in the receiver FOV, receipt of laser {wises from Earth is ensured provided: (1) 
the Earth ground station pointing error is less than the ground laser beam divergence (nominally about 50 
prad); (2) the range receiver FOV is larger than the Earth disk and is adequately boresighted with the CCD 
array; and (3) there is sufficient signal to surpass the detection threshold of the receiver. Angular errors due 
to uncertainties in planetary and most other ephemerides to important bodies in the inner Solar System are 
typically far smaller (<0.01 prad) than the nominal laser beam divergence assumed here (50 prad) as are 
the pointing control errors in a star-calibrated, meter-class telescope/tracking system (< 1 5 prad). 
Furthermore, arriving laser photons can be detected by a quadrant ranging detector, which is co-aligned 
with the center of the CCD array. As in NASA's developmental SLR2000 satellite laser ranging station, the 
quadrant detector permits fine pointing corrections of the transponder receiver at the subarcsecond level 
[Degnan and McGarry, 1997], 


4.2 TRANSMITTER POINT-AHEAD 

In a general transponder experiment, the CCD Earth image tells us where the Earth was situated one transit 
time earlier whereas the narrow beam transmitter must be pointed to where the Earth will be one transit 
time later. Thus, one transit time later and in the image plane of our receiver CCD, the Earth will he 
somewhere on the perimeter of a circle centered on our Earth image and having an angular radius equal to 
the differential pointing angle. The projection of the Earth’s forward velocity vector in the receiver plane is 
then needed to determine the precise point on the circle to which the laser beam should be directed. The 
forward velocity vector can be inferred from simultaneous CCD observations of the Moon and/or stars 
within the camera FOV. Although the Moon is significantly smaller than the Earth and therefore reflects 
fewer photons into the CCD array, it also deposits them into fewer pixels so that the contrast is not terribly 
different. This will be demonstrated in Section 5. The Moon is expected to be visible about 80% of the time 
[Mallama, 1998] and, when it is not, background stars within the nominal 1° x 1° CCD FOV can be 
substituted as available. 
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Differential pointing of the transmitter can be most easily accomplished in a bistatic optical system with a 
pair of azimuthally driven Risley prisms in the transmit path as in Figure 5. The magnitude of the angular 
offset is determined by the relative orientations of the prisms (which remains fairly constant during a 
transponder experiment), and the direction is changed by driving them in unison to locate the precise point 
on the circle. As in Figure 5, transmitter feedback to the receiver can be provided by a prism or two mirror 
reflector which samples the beam as it exits the transmit telescope and reflects it into the receiver for time 
logging and/or verification of the transmit angle. The bistatic optical design helps to eliminate laser 
backscatter from the local atmosphere as an important noise source. 


4.3 ACQUISITION AND TRACKING OF THE OPPOSITE TERMINAL IN RANGE SPACE 

To simplify the discussion of acquisition in range space, we will assume that the two transponders have a 
common laser fire interval, r qs . During each interval or cycle, each terminal records the time of departure of 
the outgoing pulse and the time of arrival of any incoming pulses relative to a fixed sub-clock, which is 
counting at the nominal laser fire rate and is derived from a highly accurate frequency standard such as a 
rubidium or cesium clock. In the case of a Mars lander, an onboard knowledge of the Earth and Mars 
ephemerides and planetary spin axis orientation and rotation rates permits a fair a priori estimate of the 
interplanetary distance (typically within a few tens of kilometers), the differential pointing angle (within an 
arcsecond) and the forward direction of the Earth’s motion. If the distance between terminals can be known 
a priori to well within the distance light travels in one laser fixe interval, the incoming pulse can be easily 
matched up with the proper cycle at the opposite terminal where it originated Because light transit times 
within the inner Solar System can span several tens of minutes, there may be thousands to millions of 
cycles between the matching cycles for each terminal, depending on the laser fire rate and interplanetary 
distance. During initial acquisition, the initial uncertainty in the laser time of fire at the opposite terminal 
can be as large as the laser fire interval itself if the clocks at Terminals A and B have not been compared in 
some time and are totally uncorrelated. Once successful two-way transponding is established, however, 
this uncertainty shrinks rapidly, allowing the receivers to be gated over a much narrower time interval. 

Every range receiver is characterized by a range gate , a range bin , and a detection threshold. In a single- 
ended range measurement to a passive reflector, the range gate , provides a temporal filter for the 
reduction of background noise and is chosen large enough to encompass the uncertainty in our a priori 
knowledge of the range to the target. The range bin y is the time over which the receiver integrates the 
incoming signal and is optimally chosen, in conventional high SNR systems, to be just large enough to 
capture the majority of the signal photons within a single pulse. In the case of photon-counting receivers 
[Degnan, 2000b], a somewhat larger range bin is chosen to collect the photons from a collection of multiple 
pulses over a time interval referred to as a frame. In either case, the number of range bins within the range 
gate is given by the simple formula 


We refer to the 2D areas defined by the horizontal borders of the “range bins” and the vertical borders of 
the “frames” in Figure 7 as cells. A consecutive sequence of frames is a super-frame [Degnan, 2000b], We 
tentatively identify potential signal cells within a given frame by counting all of the photoelectrons 
generated within each cell and comparing it to a frame threshold. We refer to this basic data processing 
scheme, which was used in LLR data processing [Abbott, et al, 1973], as Post-Detection Poisson Filtering 

In a satellite laser ranging (SLR) system, the range between ground station and satellite changes rapidly and 
the position of the range gate is varied in real time, based on an a priori range estimate, in an attempt to 
keep the satellite return centered in the gate. The measured pulse times of flight, displayed on the ordinate 
axis relative to the center of the range gate versus elapsed time on the abscissa, is referred to as an 
Observed Minus Calculated (O-C) plot [Degnan, 1985], If our a priori range and range rate corrections and 
satellite force models were perfect and additionally there were no range or time biases in the computed 
orbit, the corresponding “Observed Minus Calculated” (O-C) curve would place all of the observed signal 
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photons within a single range bin centered in the range gate. In this highly idealized example, the temporal 
width of the signal data distribution would then be determined by the timing precision of the range receiver 
and small atmosphere-induced fluctuations which, in modem SLR systems, is characterized by a one sigma 
RMS single-shot range scatter of one cm (67 picoseconds) or less [Degnan, 1993]. Background photons 
and detector dark counts, on the other hand, are randomly distributed throughout the entire range gate. 
Thus, when histogrammed according to the individual range bins, the signal counts stand out from the 
background due to their high "temporal” correlation in the range axis [Degnan, 2000a]. 

In real SLR systems, range biases displace the signal from the center of the range window, and orbital time 
biases introduce a slope in the signal data as viewed in an OC plot. Once a slope is observed in the data, 
however, a time bias correction can be iteratively applied in the range model as necessary to reduce the 
slope and ultimately the size of the range bin and/or gate for better signal contrast and noise rejection. A 
similar phenomenon will occur in an asynchronous transponder link when a range-rate, estimated from 
planetary ephemerides and applied to the photon arrival times, has a residual error. As previously 
mentioned in Section 2. 1, the error in range rate computed from contemporary ephemerides for Earth and 
Mars is estimated to be less than 1 cm/sec [Lemoine, 2000] corresponding to a slope of less than 33 
psec/sec in the O-C curve. The component of range rate due to the rotation of station A about its planetary 
axis is a sinusoidal function with peak amplitude 


u 


= m A r A cosX A cos i A 


( 39 ) 


where tu A is the angular rotation rate about the axis, r A is the volumetric mean radius of Planet A, X A is the 
latitude of the station, and i A is the inclination of the spin axis to the orbital plane. A similar expression 
applies to the station on Planet B. Using the constants in Table 1, the maximum range rate introduced by 
the rotation of Earth and Mars for equatorial stations is 67.6 m/sec and 34.6 m/sec and decreases with 
latitude. Thus, equation (39) results in OC slopes between 0 and 230 nsec/sec due to Earth rotation and 
between 0 and 120 nsec/sec due to Mars rotation. The maximum values are large enough to suggest that an 
a priori correction for planetary rotation be applied to the data prior to producing the OC plot, especially if 
frame times are longer than about 0.01 sec. 

Double-ended transponder systems use two clocks and two lasers whereas a single-ended SLR system uses 
only one laser and one clock in recording the pulse time of flight, and this difference leads to some 
additional range errors as recorded in a transponder OC plot. For example, if the frequency standards 
(clocks) at the two interacting transponder terminals have a frequency offset, there will be an additional 
contribution to the observed slope given by the familiar Doppler equation, i.e. 



( 40 ) 


where/ is the nominal clock frequency at both terminals and Af is the offset frequency between the clocks 
as discussed in Section 2.3. A fractional clock offset of 1 x 10* n corresponds to a range rate error of less 
than 3 mm/sec, or about 20 psec/sec, and is comparable in magnitude to the a priori range-rate error from 
the ephemerides. This level of offset or smaller is easily achieved by a rubidium or cesium standard, but, 
for an undisciplined quartz crystal oscillator, the slope could be up to two orders of magnitude higher. 

A further difference arises from the use of two lasers. In a single-ended SLR system, the RMS scatter of the 
data depends on the jitter in the detector response, the resolution of the timer, and the impulse response of 
the satellite, but the jitter in the laser fire time cancels out of the range measurements since the start and 
stop pulses input to the interval timer are affected equally. However, in a two-way asynchronous 
transponder, the unequal jitter in the fire times of the two lasers leads to an additional random range error, 
which is expected to dominate the overall jitter contribution and broaden the signal data in the O-C plot 
relative to die single laser case. 
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Any slope or breadth to the signal data must be accommodated by our choice of range bin. The breadth or 
precision of the signal data sets a lower limit on the range bin size, even for automated high SNR single 
pulse detection. With a residual data slope, any increase in the receiver integration (frame) time in a low 
SNR photon-counting system must be accompanied by a proportional increase in the range bin width to 
ensure that we capture all of the signal photons from multiple laser fires in a single cell. One rule of thumb 
that appears to work well in simulations is to choose a bin size at least twice as large as the expected overall 
RMS variation in the O-C range data over a frame, i.e. 

h * 2 ^(( cr e 2 ) + ( °v 2 ) + ( ~°c )) + (d ) + (^ ) (41) 

where <a 2 >, <a 2 > y and <cr c 2 > are the expected variances in the residual slopes in O-C space due to errors 
in ephemerides, planetaiy rotation rates, and clock offsets respectively, tf is the frame time, and <Xa> and 
<t b 2 > are the variances in the recorded fire times (including laser, detector and timer jitter) at terminals A 
and B respectively. 

The detection threshold, K, is used to determine the probable presence of signal against a noise 
background For single pulse detection in a high SNR system, the threshold is usually set by hardware 
whereas, for low SNR photon-counting systems, the threshold is usually set by a combination of hardware 
(e.g. a range receiver or multichannel scalar) and software which compares the counts in each bin/cell to a 
frame threshold, K, and rejects counts below the threshold as probable noise [Degnan, 2000b]. In either 
case, choosing too high a threshold results in the loss of valid range returns, whereas choosing too low a 
threshold results in increased noise-induced false alarms. Furthermore, because of the potentially large 
number of cells in a frame (especially during acquisition when range uncertainties are largest), it is possible 
that a sizable number of noise cells in a given frame will be falsely identified as signal even when the 
probability of false alarm for any given noise cell is relatively small. One approach to optimizing the frame 
threshold is to maximize the Differential Cell Count [Degnan, 2000b], defined as the mean number of 
correctly identified signal cells minus the mean number of false alarms in cells containing only solar 
background counts within a superframe consisting of M frames, i.e. 

SN=M[p;-N t ,,F’'\ <«) 

where Pf is the probability of correctly identifying the signal cell within a frame and P/ a B is the probability 
of falsely identifying a noise cell as signal. This approach yields an "optimum" threshold condition for 
Terminal B given by [Degnan, 2000b] 


K B Nj+MjLnl - (Cs zML + ln(^,J (43) 

opt In C B In C B 

where N bm is the number of range bins defined by (38), N, B is the mean signal count in the signal cell, N b B is 
the mean noise count in any given cell, and Cg is the signal cell contrast defined in the next subsection. 
From Poisson statistics, the probability of correctly detecting the signal cell is given by 
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where, in the limit of "large” mean counts (>15), one can use the Central Limit Theorem to approximate the 
Poisson distribution by a Gaussian (normal) distribution and erf(x) is the familiar error function. Similarly, 
the probability of false alarm is given by 
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The error function expressions can be used in the usual case where the optimum threshold lies between the 
peaks of the noise and signal cell populations, i.e. if we satisfy the condition N b B < < C B N b B = N 3 B 

+N b B 

For weaker signals, a powerful second test can be applied, if necessary, by requiring that cells in adjoining 
frames be correlated. In the most general sense, this simply means that cells tentatively identified as 
containing ‘‘signal” in adjacent frames must obey applicable physical laws or constraints. For example, the 
physical laws governing planetary or spacecraft motion do not allow the transponder position to make 
unexpected discontinuous jumps into widely separated range bins between frames. Thus, we can define a 
“valid trajectory” as one where the planet or spacecraft position changes by no more than one range bin in 
moving between frames and monotonically moves in the correct direction on either side of the range 
extrema. This correlation requirement allows us to apply an “N of M” test on multiple cells which survive 
the initial threshold test within a superframe and recover missing signal cells . For signal verification, the 
“N of M” test requires that at least N cells, all satisfying the threshold criteria and lying on a valid 
trajectory, be detected within the M successive frames comprising the superframe. Based on successful 
application of this test, any signal data in the up to (M-N) signal cells, which may have originally failed to 
meet the threshold criteria, can be successfully restored via interpolation between frames [Titterton et al, 
1998; Degnan, 2000b]. 

To summarize, the search in the third dimension, range, is carried out by: (1) breaking up the laser fire 
interval (or some smaller range window depending on the level of our a priori knowledge) into 
appropriately sized range bins; (2) applying corrections to the raw photon arrival times using a priori 
information on range rate to form an OC plot; (3) binning and histogramming these corrected ranges; (4) 
applying a threshold test to each bin/ceil to determine candidate signal cells; and (5) ultimately applying a 
software filter to make a final selection of signal cells, calculate and correct for any residual range-rate 
error to further improve the signal contrast and narrow the range bin/window for better noise rejection. In 
prior publications [Degnan and McGarry, 1997; Degnan, 2000a; Degnan, 2000b], we have referred to the 
combination of hardware and software that performs these functions as a Correlation Range Receiver 
(CRR). Post-Detection Poisson Filters and Correlation Range Receivers, as they apply to laser altimeters 
ranging to a planetary surface from orbit, are discussed in depth in a companion article in this journal 
[Degnan, 2000b]. Such filters tend to be simpler and even more effective when used in single-ended 
ranging or double-ended transponder systems because, as mentioned previously, the signal data falls along 
a straight line in OC space whereas surface topography results in a confined but random path over short 
time intervals. 

4.4 SIGNAL CONTRAST AND SCALE ENERGY 

The signal cell contrast , C Bi at Terminal B is defined as 


c b - 1 + 
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and n, B is the mean number of received signal photoelectrons per laser fire, nf is the mean noise 
photoelectrons per range bin, and E s f is the scale energy for Terminal A. Using (7) and (8) in (46), we can 
derive the following expression for the scale energy 
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where the sum is over the relevant noise count rates described in previous subsections. We also note from 
(47) that there is a scale energy, /?**, for each noise source, and they are additive. 

In the limiting case of single pulse detection where rt is approximately equal to the full laser pulsewidth, 
the latter equation yields a scale peak power for transmitter A. 



(48) 


which yields a signal contrast of C B = 2 at Terminal B. 

Let us first consider a transponder which is either enroute to or in orbit about a target planet or is on the 
surface of a planet or other body having no atmosphere. The scale energy for the Earth terminal (A) is then 
dominated by planetary albedo, which the outgoing laser signal from Earth must overcome, plus any 
instrumentation noise, such as detector dark counts. From (27) and (47), the scale energy for solar 
scattering from the planetary surface and atmosphere (albedo) is given by 
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which is the Transmitter A pulse energy which produces the same number of signal counts in the signal cell 
as the background noise, thereby yielding a signal cell contrast of 2 at Terminal B. 

For a transponder on the surface of the planet, we must allow for additional noise sources. We can use (19) 
and (47) to compute the scale energy for local atmospheric scattering of solar (or lunar) radiation, obtaining 
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where the approximation holds if the atmosphere of Planet B is highly transmissive. In this instance, the 
source of the exoatmospheric spectral irradiance, Nx, is the Sun during local daylight transponder 
operations or possibly one or more moons during local night operations. Under most practical situations, 
the noise contribution of lunar scattering from the atmosphere tends to be negligible when compared to the 
planetary albedo. 

Finally, we can compute the time-dependent scale energy for laser backscatter in the local atmosphere from 
(34) and (47), i.e. 
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It should be noted that the aforementioned scale energies for the three noise sources depend on only four 
instrument parameters - the transmitter A solid angle, , the receiver B solid angle, C2 r ,; the bandwidth of 
Terminal B’s spectral filter, zU fl ; and the range bin, r b . For Terminal B, the corresponding scale energies arc 
obtained by simply interchanging A and B in the above equations. 
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Clearly, noise reduction is an important first step in effective transponder design. Specifically, one must : 
(1) limit the spectral width of the bandpass filter ( spectral filtering ); (2) reduce the receiver field-of-vicw 
(FOV) to the minimum value which encompasses the transmitter beam divergence and its potential 
pointing variations and/or coalignment errors with the receiver (spatial filtering), (3) gate the receiver so 
that only photons within a given time window about the expected signal return are accepted (temporal 
filtering ); and (4) choose an optimum detection threshold (amplitude filtering). In low SNR photon 
counting systems, amplitude filtering cannot be used except in a statistical sense over multiple laser fires 
and must therefore be replaced by a Post-Detection Poisson Filter (PDPF) [Degnan and McGarry, 1997, 
Degnan, 2000b]. 

5. EARTH-MARS TRANSPONDER LINK ANALYSIS 

5.1 COAXIAL, COPLAN AR, CIRCULAR ORBIT APPROXIMATION FOR AN EARTH-MARS TRANSPONDER LINK 

The expressions derived in this section are valid for more general planetary ephemerides. However, in 
generating our numerical results, we will ignore the small orbital eccentricities (0.0935 for Mars and 
0.0167 for Earth) and small relative tilt (-1.85 °) in the planetary orbital planes and assume coplanar (xy 
plane), circular, and coaxial orbits of Planets A and B about the Sun having radii R A and R B respectively. 
We also adopt the convention R B >R A .ln this approximation, the vector positions of Planet A and B at any 
given time are given by the simple formulae 
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where co A and a> B days are the orbital angular frequencies for the Earth and Mars respectively, and the time 
origin is arbitrarily chosen when Planets A and B are at the Point of Closest Approach (PCA) on the same 
side of the Sun. The interplanetary range, range-rate, and range-acceleration vectors are then given by 
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respectively. The latter vectors have magnitudes equal to 
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and 
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is the synodic angular frequency and t s is the synodic period which, for Earth and Mars, corresponds to 
approximately 780 days. 

5.2 Doppler shift and spectral filter transmission 

The range rate term introduces a variable Doppler shift in the laser wavelength as viewed by the opposite 
terminal, given by 


Ass umin g a Lorentzian shape for the transmission function of the spectral filter yields a reduced receiver 
transmission at Terminal A given by 
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where AX a is the FWHM bandwidth of the spectral filter. 

5.3 Transmitter Point-Ahead Angles 

The point-ahead angle for the transmitter at Terminal A is given by the equation 
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where tf and t* are the "forward" and "backward" one-way flight times of a photon leaving and arriving at 
Earth at a time t respectively. In terms of general ephemerides, these photon flight times are the solutions to 
the equations 
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and 
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where, in computing the pulse time of flight, we have included both the planetary Center of Mass motions 
about the Sun as well as the terminal motions about the planetary axis of rotation as discussed in Section 
4.3. The latter equations can be solved by substituting the instantaneous range at time t divided by the 
speed of light into the RHS of (59) and continuing the iteration until the desired accuracy is achieved. 

For the purposes of calculating transmitter point-ahead, we can ignore planetary rotation and use the first 
order vector Taylor expansion 


R B (t + At) = R B {t) + AtR B {t), 


(60) 


to obtain approximate expressions for the forward and backward photon times of flight between Earth and 
Mars, i.e. 
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modifies the calculated photon flight time between Earth and Mars by less than one part in 10 4 . Equations 
(61a) through (61c) can now be used in (58) to compute the transmitter point ahead angle at Terminal A. 
As always, the transmitter point ahead angle for Terminal B is obtained by interchanging the A and B 
superscripts/subscripts in the latter equations. 

5.5 ASPECT ANGLES 

The aspect angle of the solar-illuminated surface of Planet B (Mars), as viewed from Planet A (Earth), is 
given by 
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where R(t) = R B ( t) - R A {t) . Similarly, the aspect angle of the solar-ifluminated surface of Planet A, as 
viewed from Planet B, is 
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In Figure 8, we plot: (a) the Earth-Mars range, (b) range-rate, (c) range acceleration, (d) Doppler shift, (e) 
point-ahead angles, and (f) the aspect angles as a function of the synodic phase (day) for the Earth-Mars 
transponder link. Table 1 summarizes the planetary parameters used in the calculations. The maximum 
range-rate of approximately 55 km/sec and range acceleration of 8.2 mm/sec 2 both occur midway in the 
synodic period at the point of maximum range (2.52 AU). The Doppler-induced wavelength shift is 
confined to a range of about + 0.25 angstrom, which sets a lower limit on the bandwidth of a fixed spectral 
filter although tunable filters are possible. From Figure 8e, the transmitter point-ahead for the Earth station 
varies between about 25 and 33 arcseconds and is a maximum at the points of minimum and maximum 
range. For the Mars station, the curve has the same peaks and valleys, but the point-ahead angle vanes over 
a much wider range, between 0 and 4 1 arcseconds. From Figure 8f, we see that Mars is always at least 
partially illuminated when viewed from Earth while Earth, as viewed from Mars, runs the entire gamut 
from being backlit at the point of minimum interplanetary range to fully illuminated at the point of 
maximum range. 


PLANETARY PARAMETER 

EARTH (A) 

MARS(B) 

MOON (A’) 

Mean Distance from Sun, R (AU ~ 150 x 10 6 km) 

1.0 

1.52 

1.0 


365.256 

686.98 

NA 

Length of day, hours 

24 

24.657 

NA 

Obliquity of Spin Axis to Orbital plane, deg 

23.45 

25.19 

NA 

Mean Volumetric Radius, km 

6371 

3390 

1738 

Mean Surface Reflectivity (8) 532 ran, p 

0.15 (est.) 

0.15 (est) 

0.12 

Atmospheric Transmission @532 nm, T 

0.8 

0.9 (est) 

1.0 


Table 1: Planetary parameters assumed in the link calculations. 

5.6 ANGULAR ACQUISITION BY THE CCD ARRAY 

We assume in this section that acquisition of planetary and lunar images against a dark sky is relatively 
easy and that we only need to concern ourselves with daylight acquisition. The contrast of the planetary 
images against the solar-illuminated atmosphere during daylight operations on Planet B’s surface is given 
by (37) with the corresponding equation for the opposite terminal given by an interchange of A and B in 
that expression. As mentioned in Section 4. 1, since the Earth and Moon have the same solar aspect angle 
when viewed from Mars, we can use the same expression to compute the lunar contrast (excluding rare 
eclipse events). 

We now want to integrate the CCD scene long enough so that we have a high probability of correctly 
identifying the planetary image. We can use a variation of the Differential Cell Count scheme described in 
Section 4.3 to derive an optimum pixel threshold given by 

(C^-lJiV'+ln 

K ™ = 

where N pa is the total number of pixels in the CCD array, N A is the number of pixels illuminated by the 
planetary image as given by (35), and 
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is the mean noise count per pixel due to solar scattering off the atmosphere collected during an integration 
time, r/. In deriving (64), we have used (19) and allowed the CCD quantum efficiency , the optical 
throughput field of view D cc d, and spectral bandpass AA cc f of the CCD channel to be different from 
the range receiver channel. Since the CCD spectral bandpass filter acts equally on the planetary albedo 
"signal" and the solar scattering "noise", the contrast of the planetary image (unlike the ranging signal) is 
unaffected by our choice of spectral bandwidth. Furthermore, the effective differential pixel count will 
asymptotically approach unity as the mean noise cell count increases, and we should therefore choose a 
large spectral bandwidth filter and the longest CCD integration time consistent with the loop bandwidth 
required to keep the planetary image centered within the CCD array. 

5.7 Scale Energies and Power for the Earth-Mars Link 

We can use (49) and (50) to compute the scale energies for the planetary albedo and solar scattering 
contributions to the noise as a function of synodic phase. The results are plotted in Figure 9 for (a) the 
Earth transmitter (blue) and the Mars transmitter (red). We assume a range bin of 1 nsec (single 
subnanosecond pulse detection) and a zenith atmospheric transmission at Mars of 0.9. The scale energy 
scales proportionally with the range bin width. In all cases, the scale energy for both terminals (the solid 
curves in Figure 9) is always less than 0. 1 jJL Not surprisingly, the scale energy for the Mars transmitter is 
completely dominated by the need to overcome heavy solar scattering in the Earth’s atmosphere. For the 
Earth transmitter, Earth albedo dominates local solar scattering in the Martian atmosphere except near the 
middle of the synodic cycle when the Earth is at its farthest distance from Mars. 

5.8 high SNR example: NASA MOBLAS station ranging to a "Green MOLA" 

Our first link example will be a conventional high Signal-to-Noise Ratio (SNR) link between a NASA 
Mob ile Laser (MOBLAS) station [Degnan, 1985]and a 532 nm (green), subnanosecond version of the Mars 
Qrbiter Laser Altimeter (MOLA) instrument [Ramos-Izquierdo et al, 1994], both of which are operating at 
a zenith angle of 30°. Figure 10 shows a photo of the two terminals. The laser energy of the "green MOLA" 
was chosen to provide a balanced transponder link with MOBLAS, as discussed in Section 2.5, assuming 
the 50 cm MOLA receive aperture. Coincidentally, the green pulse energy of 43 mJ, multiplied by the 
MOBLAS repetition rate of 5 Hz yields an average power of 215 mW which is roughly half the 1064 nm 
power produced by the MOLA instrument at 10 Hz. Thus, the reduction in repetition rate compensates for 
the roughly 50% loss in laser power due to conversion of the fundamental Nd:YAG wavelength from 1064 
nm to 532 nm so that the "green MOLA" should consume roughly the same prime power as the original 
MOLA. In both systems, we assume an improved transmitter beam divergence of 50 jirad with a somewhat 
wider receiver FOV of 100 jorad in order to accommodate residual boresight errors/drifts between the 
transmitter and receiver. Table 2 summarizes the assumed instrumental parameters for the high SNR link. 

In this high SNR example, we assume that each terminal detects single subnanosecond pulses, i.e. one laser 
fire per cell. We therefore select a range bin size, r b = 1 nsec, which comfortably accomodates all of the 
signal photons in a single pulse. We further assume the worst case scenario of totally unsynchronized 
clocks at Terminals A and B so that the receiver is initially ungated for the entire interpulse period of 200 
msec. This results in a very large number of range bins, i.e. N bin = 2x 10 8 . Using (42) through (45), we can 
plot, in Figure 11, (a) the optimum threshold; (b) the differential cell count; (c) the probability of detecting 
the signal; and (d) the mean false alarms per laser fire for the Earth and Mars terminals as a function of the 
synodic phase for both night and day operations. 
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Parameter 

MOBLAS (A) 

"Green MOLA" (B) 

Transmitted Pulse Energy, mJ 

100 

43 

Repetition Rate, Hz 

5 

5 

Average Laser Power, mW 

500 

215 

FWHM Pulsewidth, psec 

150 

<1000 

FWHM Beam Divergence, j^rad 

50 

50 

Telescope Diameter, cm 

76 

50 1 

Detector Quantum Efficiency. % 

12% 

12% 

Receiver Optical Efficiency.% 

40 

40 

Receiver FOV, jarad 

100 

100 

Spectral Bandwidth, nm 

0.3 

0.3 

Range Gate, msec (ungated) 

200 

200 

Range Bin, nsec 

1 

1 

Spacecraft Zenith Angle, deg 

30 

30 


Table 2: Instrumental parameters assumed in the link calculations for a conventional high SNR 
balanced transponder pair operating between Earth and Mars. 

During the first quarter of the synodic cycle, the optimum threshold drops from roughly 12-15 pe to 2-3 pe 
then levels off for two quarter cycles before climbing again to the 12-15 pe level during the last quarter. 
Near minimum interplanetary range, the signal is higher and the optimum threshold tends to higher values 
so that the differential cell count algorithm can suppress the numter of false alarms to a negligible number 
without significantly affecting the signal detection probability. Near maximum range, the algorithm lowers 
the threshold to allow a higher probability of detection for the weaker signal but in the process makes false 
alarms more likely. The probability of successfully detecting the signal is essentially unity for the first and 
last quarter cycle for both terminals under both day and night conditions. A negligible number of false 
alarms occur (<1 per fire) under both day and night scenarios during the first and last quarters. Near the 
point of maximum interplanetary range, however, the detection probability for daylight operations drops to 
about 94.5% for the Earth terminal and 98.5% for the Mars terminal. Furthermore, because of the lowered 
threshold near the point of maximum range and the huge number of 1 nsec range bins within the ungated 
laser fire period of 200 msec, the mean false alarms per fire can jump to 10 at Earth and to abouHOO at 
Mars in spite of an extremely low false alarm probability per bin which is on the order of 5 x 10 for the 
Earth terminal and 5 x 10' 7 for the Mars terminal. Nevertheless, the false alarms are widely scattered 
throughout the laser fire interval whereas the 5 Hz signal counts will stand out prominently against the false 
counts when displayed in a histogrammed O-C plot because of their temporal coherence. Furthermore, once 
the signal is acquired, the range gate can be narrowed sufficiently to suppress virtually all subsequent false 
alarms. 

5.9 LOW snr example: SLR2000 ranging to a "green mla" 

We now provide an example of a low SNR link where multiple single photon returns in a high repetition 
rate link are accumulated in the cells of a correlation range receiver. The potential advantages of such a 
system were described in Section 2.6. The Earth ground station is assumed to be NASA's developmental 
SLR2000 station [Degnan and McGarry, 1997] and the Mars terminal is assumed to be similar in size to the 
Messenger Mercury Laser Altimeter (MLA), which is smaller and more lightweight than the Mars MOLA 
instrument. Both instruments are shown in Figure 12 and are assumed to use compact, low energy, 
microlaser transmitters which operate at a much higher 2 kHz rate [Degnan and Zayhowski, 1998], The 
laser energy of the "green MLA" was again chosen to provide a balanced transponder link with SLR2000 
assuming the 25 cm MLA receive aperture. The instrument characteristics arc summarized in Table 3. 

In this low SNR example, we select a range bin size, Tb = 2 nsec, which comfortably accomodates all of the 
signal photons. We again assume the worst case scenario of totally unsynchronized clocks at Terminals A 
and B so that the receiver is initially ungated for the entire interpulse period of 0.5 msec, but the higher 
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repetition rate results in a much smaller number of range bins, i.e. I N bm - 2.5 x 10 s . We plot, in Figure 13, 
(a) the optimum threshold; (b) the differential cell count; (c) the probability of detecting the signal; (d) the 
mean false alarms per laser fire; and (e) the number of interplanetary range measurements per second for 
the Earth and Mars terminals as a function of the synodic day for both night and day operations. Because 
we have assumed a ’"balanced" system, the frequency of recorded ranges is the same at both terminals and 
the two curves in Figure 13e overlap. 

The behavior of the high and low SNR links as a function of synodic phase are both qualitatively and 
quantitatively different. Instead of decreasing during the first quarter of the synodic cycle and leveling off 
as in the high SNR case, the optimum threshold in the photon-counting example climbs, under local 
daylight conditions, from roughly 10 pe to perhaps 30 pe for the Mars terminal and to 100 pe for the Earth 
terminal at maximum interplanetary range and then falls off symmetrically during the second half of the 
synodic cycle. This behavior occurs because, at longer interplanetary ranges, it takes longer to collect the 
signal photons necessary to unambiguously distinguish the signal cell from the noise cells. The longer 
frame time results in a higher mean noise count per cell, which in turn drives the optimum frame threshold 
higher. The effectiveness of our optimum threshold algorithm is demonstrated by the fact that the 
Normalized Differential Cell Count, with a maximum value of 1, never drops below 0.96 at either terminal 
during either day or night operations. Similarly, the probability of correctly identifying the signal cell is is 
greater than 98.7% under all conditions, and the mean false alarm rate is less than one cell in 30 frames 
under the worst case condition of unsynchronized clocks. As in the high SNR case, the number of false 
alarms is driven down dramatically once the opposite terminal is acquired and the range gate shrunk 
accordingly. We also note that the rate of range returns at maximum range is 4 Hz, only slightly lower than 
the nearly fixed 5 Hz rate of the high SNR case, but rises to approximately 100 Hz near minimum range. It 
should be further noted that the Mixed Power-Aperture Product (defined in Section 2.5) for this low SNR 
example is 0.013 Watt-m 2 or 7.5 times smaller than the same product in the high SNR case (.098 Watt-m 2 ). 
Thus, as postulated in Section 2.6, a significantly higher range sampling rate can be obtained with the low 
SNR photon-counting system over the full synodic period. 


Parameter 


"Green MLA" (B) 

Transmitted Pulse Energy, mJ 

0.130 

0.053 

Repetition Rate, Hz 

2000 

5 

Average Laser Power, mW 

260 

103 

FWHM Pulsewidth, psec 

200 

<1000 

FWHM Beam Divergence, firad 

50 

50 

Telescope Diameter, cm 

40 

25 

Detector Quantum Efficiency, % 

12% 

12% 

Receiver Optical Efficiency, % 

40 

40 

Receiver FOV, pxad 

100 

100 

Spectral Bandwidth, nm 

0.3 

0.3 

Range Gate, msec (ungated) 

0.5 

0.5 

Range Bin, nsec 

2 

2 

Spacecraft Zenith Angle, deg 

30 

30 


Table 3: Instrumental parameters assumed in the link calculations for a balanced, low SNR (photon- 
counting) transponder pair operating between Earth and Mars. 

6. SUMMARY 

Interplanetary laser transponders, operating anywhere within the inner Solar System, are well within the 
current state-of-the-art. Analysis has shown that a NASA MOBLAS SLR station, operating at 5 Hz in 
conjunction with a MOLA-sized instrument at Mars, can range with near 100% probability of detection. 
Alternatively, a compact, low power Asynchronous Microlaser Transponder, comparable in size to the 
Messenger Mercury Laser Altimeter (ML A) and working in single photon detection mode in conjunction 
with NASA's developmental SLR2000 Satellite Laser Ranging station between Mars and Earth, is capable 
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of recording between 4 and 100 two-way range measurements per second. These measurements allow 
precise determination of the range and the time offset between the ground and spaceborne docks. 

Because interplanetary light travel times are several orders of magnitude longer than is typical for artificial 
satellites, the absolute accuracy of the range and clock offset measurements is dominated by the frequency 
accuracy and stability of the ground and spaceborne clocks. Errors in the range vernier or atmospheric 
propagation delays are typically at the sub-centimeter level as for artificial satellite ranging. Decimeter 
accuracy interplanetary range measurements and subnanosecond time transfer would appear to be easily 
achievable with the conceptual system described here provided the spacecraft carried a rubidium or cesium 
frequency standard. If a ground-based maser were used to govern the Earth station timing and to 
“discipline” the inferior onboard clock significantly greater accuracies might be achieved. Based on past 
experience with the Viking lander on Mars, current microwave system precisions appear to be limited by 
hardware at the few meter level, and absolute accuracies may be much worse due to large uncertainties in 
propagation delays introduced by the interplanetary solar plasma at the lower microwave frequencies. 

Potential science benefits of laser transponders would include more precise ephemerides of selected planets 
and asteroids, improved knowledge of planetary librations and their interior makeup, the ability to perform 
more precise general relativistic experiments on an interplanetary scale, the study of solar gravity and 
internal mass distribution using Sun-orbiting satellites (much as we currently measure the Earth's gravity 
field using artificial satellites) , and measurement of the mass distribution within the asteroid belt. The 
latter could be accomplished by landing a transponder on an asteroid and measuring the range variations as 
it circulates through the belt. Besides providing relief to an already overworked Deep Space Network 
(DSN) and an independent assessment of range and range-rate, engineering benefits of transponders might 
improve the reliability of interplanetary spacecraft operations through more accurate spacecraft trajectories 
and the ability to transfer time from Earth standards with subnanosecond accuracy. 

Successful field tests of the SLR2000 system in the 2001-2002 time frame, combined with our ongoing 
laboratory transponder tests and continued launches of laser-based instrumentation such as the Geoscience 
Laser Altimer System (GLAS) and Vegetation Canopy Lidar (VCL), should remove most, if not all, of the 
technical risk and uncertainty in proposing an interplanetary microlaser transponder mission to fly within 
the next few years. In addition, we have been funded to develop and demonstrate, by August 2001, a 10 
kHz Airborne Microlaser Altimeter under NASA’s Instrument Incubator Program [Degnan, 2000b]. The 
combined technical and scientific results from these varied programs should greatly enhance the credibility 
and scientific desirability of dual-purpose instrumentation like the combined Microlaser Altimeter 
Transponder (MAT) described in Figure 5. 
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Figure 1: Timing diagram for an asynchronous laser transponder on Mars exchanging pulses with an Earth 
station. A pulse leaves the Earth station at time tEi, travels to Mars in a time tgM, and arrives at Mars at 
time tM2- Similarly, a pulse is transmitted from Mars at time t^i, travels to Earth in a time t^, and arrives 
on Earth at time t^. Times of departure and arrival are referenced to that terminal's clock, and the Mars 
values are transmitted to Earth over a microwave or laser communications link. The quantity t is the actual 
time interval between the departure of pulses from the Earth and Mars stations and, combined with the 
outgoing pulse times as recorded by their respective clocks, provides the offset between the ground and 
spacecraft clocks. 



Figure 2: Range sampling rate as a function of the detection threshold, n t , and the mean signal count, n*, 
normalized for a particular Mixed Power- Aperture Product and range. 





Terminal B Zenith Angle (deg) 

Figure 3: Normalized contribution of local atmospheric scattering to the noise at Terminal A as a function 
of the Solar and Terminal B zenith angles, Os (0°> 30°, 60°) and Ob (0 to 90 ), respectively as computed 
using (19) in the text. A standard clear Earth atmosphere (visibility = 23 km) was assumed with a zenith 
transmission of about 70%. The solid line represents the approximation given by (21) which is independent 
of the solar zenith angle and valid under the assumption of either high atmospheric transmission (7 b ~1) or 
0s ~0 b> The approximation is fairly good for 0b < 60° and improves for clearer or more transmissive 
atmospheres. 
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Figure 4: (a) Solar system view of Planets A and B circling the Sun. The aspect angle for Planet A is 
defined as the angle formed by two vectors originating at Planet A and directed toward the other two bodies 
and vice versa, (b) The fraction of the planetary disk illuminated by the Sun as a function of aspect angle, 
(c) The fraction of the Earth and Mars disk illuminated by the Sun as viewed by the transponder terminals 
as a function of synodic phase (see Section 5). In our modeling of planetary albedo noise, each planet is 
assumed to be a solid sphere whose surface is characterized by a globally averaged surface reflectivity, p, 
with a thin atmospheric shell characterized by a zenith transmission, T . Planetary albedo noise is generated 
from any point on the planet where the solar illumination and receiver FOV overlap. 





MICROLASER ALTIMETER TRANSPONDER (MAT) BLOCK DIAGRAM 



Figure 5 : Block diagram of a dual-mode Microlaser Altimeter Transponder (MAT) instrument. The CCD 
array is used to crudely center the planetary image within the quadrant photomultiplier FOV so that the 
latter can provide subarcsecond pointing corrections to maximize the link. A Risley prism pair is used to 
offset the transmitter line of sight for "point-ahead". The accuracy of the point-ahead can be monitored by 
leaking a small portion of the outgoing pulse into the receiver and onto the CCD. The transponder can be 
used at any point with in the mission, i.e. interplanetary cruise, in orbit, or the surface of another planet. 
While still in orbit, the same system can produce accurate 3-D topographic maps of the surface prior to 
landing and record coaligned 2-D visual images using the high resolution CCD camera. 
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Figure 6: (a) Simulated CCD camera display of a 1° by 1° FOV showing the centered Earth, Moon, and 
background stars used in the computation of the transmitter point ahead angle and direction, (b) The circle 
around the Earth represents the magnitude of the angular offset as computed from the planetary 
ephemerides, and the arrow represents the Earth’s forward motion as determined by other celestial objects 
(Moon, stars) in the camera FOV. The intersection of these two curves indicates the desired laser beam 
pointing. Sampling of the outgoing laser beam by the prism assembly can provide pointing verification. 

(c) Number of CCD pixels illuminated at the Earth (blue) and Mars (red) terminals by the partial disk 
image of the opposite planet; dashed red is the Moon as observed from Mars, (d) Contrast of planetary 
image to background counts caused by solar scattering from the local atmosphere during daylight 
operations. 
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Figure 7: The arrival times of incoming photons are recorded relative to when they are expected in an 
Observed Minus Calculated (O-C) plot In this reference system, a Correlation Range Receiver (CRR) 
divides the range gate in the vertical axis into range bins and the sequence of laser fires along the 
horizontal axis into frames The bin and frame sizes are chosen to ensure high probability that the observed 
signal counts all fall within a single cell , 2-D areas defined by the vertical borders of the frame and the 
horizontal borders of the bins. The choice of bin and frame size are intertwined and must take into account 
the basic precision of the measurement, which sets a minimum bin size, as well as the time it takes to 
collect enough signal counts to adequately discriminate the signal cell from the noise cells. Various 
residual uncertainties in the range rate, due to ephemerides and planetary rotation model errors, as well as 
clock frequency offsets, which can mimic a range-rate error via the Doppler equation, can cause the signal 
photons to fall along a sloped line with the random scatter superimposed. 











Figure 8: Various plots relevant to the Earth-Mars link as a function of synodic phase/day: (a) range in 
AU; (b) range-rate in km/scc; (c) range acceleration in cm/sec 2 ; (d) Doppler shift in Angstroms; (e) 
transmitter point ahead angle in arcseconds; and (f) aspect angles. In (e) and (f), the Earth terminal is 
represented by a solid line and the Mars terminal by a dashed line. 










Figure 9: Scale energies in microjoules as a function of synodic phase for both day and night operations 
with both terminals at a zenith angle of 30°. Scale energies for the Earth and Mars transmitters are 
represented by the blue and red curves respectively. The dashed lines represents the planetary albedo 
contribution and represents night operations. The dotted line represents the additional contribution from 
solar scattering in the local atmosphere. The solid line represents the sum of both scale energies 
representing local daylight operations. 
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Figure 10: Systems making up the sample high SNR transponder link: (a) MOBLAS Satellite Laser 
Ranging Station; (b) Mars Orbiter Laser Altimeter (MOLA) instrument. 
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Figure 11: Operating parameters as a function of synodic phase/day for the balanced, high SNR Earth- 
Mars transponder link operating at 5 Hz during initial acquisition under worst case conditions of 
unsynchronized clocks: (a) optimum detection threshold; (b) optimized differential cell count, (c) 
optimized signal detection probability ; and (d) mean false alarms per laser fire. The Earth and Mars 
terminals are represented by blue and red respectively. Daylight operations are indicated by the solid lines 
and local night operations are indicated by dashed lines. 
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Figure 12: Systems making up the low SNR (photon-counting) transponder link: (a) SLR2000 Satellite 
Laser Ranging Station; (b) Messenger Mercury Laser Altimeter (MLA) instrument. 
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Figure 13: Operating parameters as a function of synodic phase/day for the balanced, low SNR (photon- 
counting) Earth-Mars transponder link: (a) optimum detection threshold; (b) optimized differential cell 
count; (c) optimized signal detection probability; (d) mean false alarms per laser fire; and (c) the number of 
interplanetary range measurements per second (same for both terminals since they are "balanced’ ). Blue 
lines indicate the Earth Terminal; red lines indicate the Mars Terminal. Solid lines represent daylight 
operations; dashed lines represent night operations. 








